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I. INTRODUCTION 

Within the span of one year, multiple independent neu- 
trino oscillation measurements by reactor and accelerator 
long-baseline experiments have established the neutrino 
6*i3 mixing angle to be ~ 9° [THS]. For neutrino oscil- 
lation experiments, the newly gained knowledge of the 
sin 2 2# 13 mixing strength brings significant clarifications 
concerning sensitivity to CP-violating effects (CPV) as 
may ensue with either the normal hierarchy (NH) or the 
inverted hierarchy (IH) for the neutrino mass eigenstates. 

The revelation of # 13 has enabled first-time determi- 
nations of "exclusion curves" for the Dirac CP phase, 5, 
for each mass hierarchy in experiments having sensitiv- 
ity to perturbations in neutrino oscillations arising from 
terrestrial matter effects: The SuperKamiokande collab- 
oration has reported x 2 versus S fits for NH and for IH 
using atmospheric neutrino data [6] ; exclusion confidence 
level (C.L.) curves are reported by the MINOS collabo- 
ration based upon v e and V e appearance at the exper- 
iment's 735-kilometer baseline [TJ. While the exclusion 
levels thus far achieved are quite modest, they serve to re- 
mind that a new era of experimental scrutiny of neutrino 
flavor oscillations is getting underway and it remains to 
be seen whether conventional three-flavor mixing phe- 
nomenology will continue to be an adequate framework. 
Observational deviations from this framework are an ex- 
citing possibility as plausible harbingers of physics be- 
yond the Standard Model. As a source of such deviations, 
the possibility that neutrinos propagating through dense 
matter may participate in effective, neutral-current like 
nonstandard interactions (NSI) has received considerable 
attention for more than a decade [8]. 

In standard three-flavor neutrino oscillation phe- 
nomenology, the Hamiltonian in flavor basis includes the 
Mikheyev-Smirnov-Wolfenstein (MSW) matter poten- 
tial. The MSW potential accounts for coherent forward 
scattering of electron-flavor neutrinos from the electrons 
of ambient matter [S]. In an NSI scenario, the Hamilto- 
nian carries additional matter potential terms analogous 
to the MSW potential which allow flavor-changing as well 
as flavor-conserving NSI scattering processes. There are 
six possible NSI amplitudes which can arise in neutrino 



propagation through matter. These include three real- 
valued, flavor-diagonal amplitudes conventionally desig- 
nated as e ee , e w , and e rT , and three flavor-changing am- 
plitudes which may carry CP-violating phases: e ep , e er , 
and e^r. The phenomenology of neutrino NSI in prop- 
agation for accelerator, atmospheric, and solar neutri- 
nos and for various beam-plus-detector(s) configurations, 
has received extensive treatment. Data from hadron and 
lepton colliders have also been utilized in NSI studies, 
for NSI couplings involving quarks or electrons can give 
rise to anomalous monojet, monophoton, and multilcp- 
ton events [TU] . Comprehensive citations to the published 
literature can be found in [H [TTHI3] . 

The available data allow upper bounds to be set on the 
magnitudes of NSI couplings. According to the analysis 
of Ref. [H], the effective NSI parameters for terrestrial 
matter are |e MM | < 0.07, |e eM | < 0.33, and \e^ T \ < 0.33 
at 90% C.L. Additionally, on the basis of consistency 
with the high-energy atmospheric data, it is proposed 
that bounds of a few percent are appropriate for \e e/i \ 
and |e MT | [T^llIS]- For |e ee |, |c eT |, and \e TT \ however the 
bounds at 90% C.L. are much weaker; Ref. [2] finds 
these to be < 4.2, < 3.0, and < 21 respectively. For neu- 
trino NSI with electrons (but not with u or d quarks), 
more stringent limits have been set for e^, , e^?, and 
based upon analysis of solar and KamLAND neu- 
trino data [16] . With respect to NSI limits derived using 
monojet plus missing energy datasets of hadron collider 
experiments, our characterization is appropriate for the 
"light mediator" regime of anomalous monojet processes 
involving NSIs [10l[T7]. Thus, at present, the e eT , e ee , 
and e TT NSI are so poorly constrained that large mat- 
ter effects, of strengths which rival or exceed that of the 
MSW matter potential, remain as viable phenomenolog- 
ical possibilities. These three couplings enter into the 
probabilities for v e appearance oscillations, and the per- 
turbations they may introduce can be searched for by 
observing and v ^ beams at long baselines. 

Concerning e eT V e , the current situation with the solar 
8 B neutrino energy spectrum is worthy of note [T81 119j . 
The low-threshold measurements carried out thus far by 
Borexino, Super-Kamiokande, and the Sudbury Neutrino 
Observatory do not exhibit an upturn with decreasing 
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E v as predicted by standard oscillations with the MSW 
matter effect. However a deviation from conventional 
oscillations in the form of a flatter spectrum is naturally 
provided by an e eT coupling strength O(10 _1 ) [T71 |2"U] . 

In this work we examine manifestations of the NSI 
matter potential e eT V e together with e ee V e and e TT V e as 
may occur in neutrino propagation through the constant- 
density terrestrial crust. We focus on v e and V e appear- 
ance oscillations and evaluate the implications of recent 
signal event counts reported by the T2K (295 km) [T1I2T] 
and MINOS (735 km) [22] long-baseline experi- 
ments. The sensitivity of conventional, terrestrial long- 
baseline experiments to the e eT V e NSI has been explored 
in previous works by other researchers. In particular, ex- 
perimental e eT V e sensitivity has been examined for the 
295-kilometer baseline of T2K [53] , for the 735-kilometer 
baselines of MINOS and OPERA [HH2E], for the 810- 
kilometer baseline of NO^A [17], and for the 1050 km 
baseline proposed for T2KK [12]. Several of these stud- 
ies make use of constraints deduced from testing NSI 
scenarios using data from atmospheric neutrino exper- 
iments [HI [29] [30] . The work reported here utilizes the 
insights from these previous studies and examines the 
most recent observations of positive v e appearance in two 
accelerator beam long-baseline experiments in light of the 
newly delineated value range allowed to #13. Our treat- 
ment is restricted to neutral current NSI processes as 
may occur with neutrino propagation in matter. Com- 
plications arising from possible NSI effects in neutrino 
production and/or detection processes in current exper- 
iments are not considered [5]. 



II. OUTLINE 

We proceed as follows: In Sec. Ill wc define the 
matter Hamiltonian to include the NSI e eT , e ee , and 
e TT , and we assume, on the basis of bounds previously 
proposed [HI [13] that the e MM , e eM , and e MT are much 
smaller and can be neglected. We then present a formal- 
ism which characterizes three-flavor neutrino oscillations 
with NSI. Specifically, we express the v e appearance am- 
plitude A{v^ — » v e ) as a sum of three terms, Ti (z=l,2,3), 
the absolute square of which gives the exact appearance 
oscillation probability for neutrinos traversing terrestrial 
matter of constant density. Our expressions are obtained 
by deriving the time evolution operator for which the 
three-flavor Hamiltonian including matter effects is the 
generator. The methodology for this approach is pre- 
sented in Ref. [31 ; a summary of the derivation with 
inclusion of the NSI considered here is given in the Ap- 
pendix below. The analytic forms serve to illuminate the 
relative contributions arising from the various NSI and 
from the CP phases 5 and 5 eT . Their compact nature is 
effective in reducing input for computation, thereby in- 
creasing algorithm speeds. The analytic forms have been 
used to check our fits to the data which are carried out 
using numerical techniques. 

In Sec. IV we summarize the observations of the MI- 



NOS and T2K experiments that we use in order to fit for 
the NSI couplings. In Sec. V we present a sequence of 
NSI fits to the data. We commence with a minimalist sce- 
nario, namely that e er is the only active NSI and that its 
coupling is real- valued (hence neglecting its phase S eT de- 
gree of freedom). We then fit for the magnitude |e eT | and 
the sum of the CP phases S + 5 eT . Finally we consider the 
realistic situation wherein complex e eT is active together 
with the flavor-diagonal NSI couplings e ee , and e TT . For 
the latter we introduce a constraining relationship which 
is based upon the behavior of disappearance oscilla- 
tions at high energies for atmospheric neutrinos. This 
allows the number of variables in the fit to be limited 
to |e er |, e ee , together with the above-mentioned sum of 
CP phases. In Sec. VI we summarize the constraints 
for |e er | which are indicated by our fits, and take note of 
near-term experimental developments which will enable 
these constraints to be improved. 



III. EXACT AMPLITUDE FOR ->■ v e 
OSCILLATIONS 

A. Three-flavor oscillations with NSI matter effects 

For neutrino propagation in vacuum, the Hamiltonian 
in the basis of three mass eigenstates (i =1, 2, 3) is 
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diag (0, a, 1) , where 



E v 



and a = 
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are, respectively, the vacuum oscillation length and the 
mass hierarchy ratio. The transformation from mass ba- 
sis to neutrino flavor basis {f = e, /1, r) is 
provided by the unitary mixing matrix 
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wherein the atmospheric and solar mixings are accounted 
for via the rotation matrices Ri(023) and Rz{d\2)- The 
Dirac CP phase 6 is included via the auxiliary matri- 
ces la = diag(l, 1, e lS ) and I_s = Then the vacuum 
Hamiltonian in flavor basis is given by the unitary trans- 
formation 
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(v) 



Vmix^o u mix i 



(4) 



and the effective wave equation for vacuum propagation 
of flavor states is 



(•5) 



To we add (in flavor basis) the MSW and NSI mat- 
ter interactions: 
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Here, V e — splGpn,, is the MSW matter interaction [5] 
where Gf is the Fermi coupling constant and n e is the 
electron density in matter. The standard MSW matter 
effect is modified by the presence of the real- valued, di- 
agonal NSI interactions e ee V e and e TT V e , and by the off- 
diagonal e eT V e interaction. Since the latter amplitude 
may carry a CP- violating phase, S eT , hereafter we desig- 
nate the magnitude |e er | and display the phase explicitly. 
It is convenient to absorb V e into the matter potential, 



A = 2£ v V e , 



and to write Eq. Q as 
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A method to solve the time evolution operator in flavor 
basis U(v\t = £,0) for propagation to baseline distance 
I in constant density matter is presented in Ref. [3"T] . 

for conventional three-flavor oscillations with flatter — 
diag (A/2£ v , 0, 0). This same approach can be used to 
obtain an exact solution for the more elaborate matter 
interactions of Eq. |8]). Following Ref. [3T], the matrix 
elements of the evolution operator U^ v \£) corresponding 

to Hq^ + -Hotter provide the various possible three- flavor 
oscillation amplitudes. The v e appearance amplitude is 
given by element ll[ 2 ^ which can be broken out as a sum 
of three terms: 



T 1 +T 2 + T 3 



(9) 



The component amplitudes Ti comprise an analytic foun- 
dation for our investigation of NSI constraints arising 
from the recent v e appearance observations by MINOS 
and T2K. 

In the Sections to follow we specify the Ti and then 
focus upon their implications. Details concerning the 
derivation of the evolution operator U^{£) which un- 
derwrites the Ti are provided in the Appendix. 



B. Specification of the T, amplitude terms 

In order to write compact expressions for the Tj, we 
define some notations. For the mixing angles we use s,,- = 
sin Oij and cy = cos 0,j ; for the atmospheric oscillation 
phase we write 



A 



4£„ ~ Al„ 



(10) 



We define scaled forms a' and a" for the hierarchy pa- 
rameter : 



a = sin 2^i2 ■ a , and 
a" = (1 - ic\ 2 ) ■ a . 



(11) 



The mixing strengths involving 6*13 are often accompa- 
nied by the factor (1 — s\ 2 a), hence we define 



sin 20i 



(1 — s 12 a ) ' sm 2^13, an d 



cos20 13 = (1 - s\ 2 a) ■ cos 20 
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(12) 



We need to refer to the elements of the full Hamiltonian 
in propagation basis as obtained after a re-phasing of 
its diagonal elements (see Sec. IV. B of Ref. [31). The 
Hamiltonian at that stage has the form 




(13) 



Its diagonal elements are the following real- valued func- 
tions: 



Y (cos 20i3 - A[l + e ee - c\ 3 e TT ] ] 



^ 4£ 



e ee - (2s2 3 - c 23Kt] + a") , 



(14) 



We designate the complex- valued off-diagonal elements 
using lower-case letters as follows: 



/ = ^- (sin20 13 + 2c 23 \e eT \e^ 6 + 5 ^Aj , 
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5 " 44 
Then we have 
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Ti = H)s 23 -^ • sin(TVA) • e"* 5 
where 

N = A£ V -N = A£ V - [|/| 2 + Q 2 ]5. 
The second term in Eq. (|9| is 

T 2 



-i)c 23 --sin(fjA)-e lGA , 



where 



and 



fj = 4l v ■ V = 44 • [|r| 2 + \b\ 



G = 4£ v ■ G. 



21 2 



(15) 
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Of the three amplitude terms in Eq. T 3 is the 
most intricate. If the NSIs were known to be small, e.g. 
|e w <| < a, then T 3 could be neglected. However large 
NSIs are a distinct possibility and so T 3 is to be retained. 
For convenience we define two complex functions Si and 
S 2 : 



Si 



rb 

^2 ' 



(23) 
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and 



S 2 = (-») 



N 



Si 
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(24) 



- iS 



T3 can then be expressed as 

T 3 = -2s 23 • {Si cos^A) + 5 2 sin(rjA)} • sin*( 

(25) 

In summary, the three amplitude terms of Eq. ([9| are 
given by Eqs. (18 1, (20), and (25). These can be coded 
as complex functions, and the — > v e oscillation prob- 
ability can be constructed as |-4(f M — >■ v e )\ 2 . 



C. Appearance probability upon neglecting T3 

The ways in which NSI matter effects introduce distor- 
tions to conventional oscillations can be discerned in part 
by examining an approximate form for the v e appearance 
probability Viy^ — ¥ v e ). Under the assumption that all 
e w ' are relatively small, we may neglect T3 and write 
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In the last term, the ratio 



sin 2 (77 A). 



(26) 



reduces to cf 3 in the limit 



that the NSI couplings go to zero. More generally, the 
oscillation probability of Eq. (26) reduces to the three 
leading terms of the exact formula of Ref. [3T] in the limit 
that the NSI interactions are turned off. As discussed in 
Ref. [3T], these same three terms are related to the well- 
known perturbative formula of Cervera et al. (Ref. |35j : 
see also (35], [37]). One manifestation of a sizable \e eT \ 
occurs within the factor |/| 2 of the first term of Eq. ( [26] ) . 
Referring to Eq. (15), one sees that the term contain- 
ing \e eT \e l ( s+ScT ) A causes the effective mixing strength 
to deviate from sin26 , 13 , giving a dependence upon E v . 
Amplitude expressions which lead to Viv^ — > v e ) of ac- 



curacy comparable to Eq. ( 26 ) have been discussed in 
previous works [HJ [27] . 



IV. v e APPEARANCE IN T2K AND MINOS 



The occurrence of electron-shower dominated events 
with rates as predicted for ia, — ¥ v e oscillations, has re- 



cently been reaffirmed by the T2K and MINOS long- 
baseline experiments. In T2K, data exposures to the ex- 
periment's low-energy, off-axis (2.5°) beam totaling 
2.56 x 10 20 protons-on-target (PoT) have been analyzed. 
Among events having reconstructed energies less than 
1250 MeV, 10 v e charged-current event candidates are 
observed, to be compared to 2.47 background events pre- 
dicted for null oscillations [21]. The ten candidate signal 
events include six v e events reported previously by T2K 
as evidence for a relatively large 6*13 mixing angle pp. 

The recent MINOS results are based on exposures 
to the NuMI low-energy beam of 10.6 x 10 20 PoT 
in neutrino-focusing mode and 3.3 x 10 20 PoT in 
antincutrino- focusing mode. It is reported that, from 
data runs with ^-focusing, 152 candidate v e events are 
observed while 128.6 events are expected for null oscil- 
lations. (For NH with sin 2 26> 13 = 0.10 and 8 CP = 0, 
161.1 events are expected.) For running with Z^-focusing 
(reversed horn-current running), 20 (v e + v e ) candidate 
events are observed, while 17.5 events are expected for 
null oscillations. (For NH with sin 2 26*13 = 0.10 and 
Sep = 0, 21.2 events are expected.) [22] , 

For the purpose of fitting to NSI scenarios, we treat 
both experiments as counting experiments in which a 
signal has been measured over and above an estimated 
background. Errors are assigned according to sample 
statistics plus allowance for systematic errors associated 
with background estimation and signal detection. For 
MINOS we allot a conservative systematic error estimate 
of 6% [2J; for T2K we allot 15% [3[2T]. For the 295-km 
baseline of T2K and for the 735-km baseline of MINOS 
as well, neutrino propagation is confined to the Earth's 
crust, for which a density of p — 2.72 g/cm 3 is assumed. 
In fitting of three-flavor neutrino oscillations with matter 
effects, we use V e = 1.1 xl0~ 13 eV = (1/1900) km" 1 [32J. 



V. ALLOWED REGIONS FOR e eT 

We proceed with fitting of three-flavor neutrino oscil- 
lations including NSI to the T2K and MINOS v e appear- 
ance data. A log-likelihood fit of three terms is used to 
compare the observed versus expected signal rates for v e 
appearance in T2K, and for v e appearance and v e ap- 
pearance in MINOS. As previously noted, we neglect the 



NSI 



e eM , and e^r on the basis of the existing up- 



per bounds, and focus on the possible role for e er which 
may be operative in conjunction with e ee and e TT . Even 
with restriction to the latter three NSI, the number of 
degrees of freedom available to an oscillation scenario 
remains rather daunting, for the CP phases S and 5 eT 
are present together with the three coupling strengths, 
and the two possibilities for the mass hierarchy must 
be considered. Additionally the fits require values to 
be specified for the atmospheric Am^and solar A?n 2 i 
mass-squared differences and for the mixing angles 6*23, 
6*13, and 6*12- For these we use the world-average values 
and 1 a error ranges obtained for the normal hierarchy 



5 



by Ref. [35]. 

In the fits, the probabilities are computed by construct- 
ing a numeric Hamiltonian in flavor basis, solving for its 
eigensystem, and using it to propagate the neutrino am- 
plitudes. The oscillation probabilities are then assembled 
and multiplied by "event densities" constructed so as to 
yield the differential event rates predicted for null oscilla- 
tions. For each experiment, integration of the oscillation- 
weighted event density over the neutrino energy range 
probed provides the number of events predicted in the 
presence of oscillations. The prediction is then compared 
to the observed number of events using the log-likelihood 
distribution given below: 



= -2 £ ln£(^ | K bs , aj) + X 2 cnalty , (27) 



where 
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In the above expressions N p and N l obs are the predicted 
and observed number of events respectively, for experi- 
mental measurements i = 1,2,3. The systematic uncer- 
tainty of is denoted by a p and is taken into account 
by minimizing the nuisance parameter £ representing a 
fractional shift in N p . Current best-fit values for sin 2 #13 
and sin 2 023 are assigned to S13 and S23, and their un- 
certainties are given by SS13 and 5s23- For all fits re- 
ported below, marginalization is carried out for sin 2 6*23 
and sin 2 6»i 3 [55] . 



A. Real-valued e CT as sole operative NSI 

As a first step, we consider a minimalist scenario in 
which e eT is the only operative NSI which we restrict to 
be real-valued, allowing it to be positive or negative but 
otherwise ignoring its phase degree of freedom. We carry 
out two sets of fits with NH and IH treated separately in 
each set. For the first set the value zero is assigned to the 
Dirac CP phase S, while for the second set the range to 
2-7T of S is marginalized over in the fits. In both sets the 
data is fitted to v^iv^) — > v e ipe) oscillations with e eT V e 
included together with the conventional MSW matter ef- 
fect. The distributions of A X 2 = X 2 ~ X 2 (best-fit) from 
the two sets of fits are shown in Figure [T] The distribu- 
tions serve as exclusion curves, with values of e er having 
A x 2 which exceed 1.0 (2.71) being excluded at 68%(90%) 
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FIG. 1. Distributions of Ax from fitting v e and v e appear- 
ance rates reported by MINOS and T2K to v^iy^ — \ v R {y e ) 
oscillations with the e eT NSI restricted to real values. All fits 
are marginalized over the allowed ranges for the 623 and #13 
mixing angles. For the fits of the dot-dash line (NH) and 
dotted line (IH) distributions, the Dirac CP phase, 8, is set 
to zero. With S marginalized in the fitting, the bounds on 
real-valued e eT become less stringent as shown by the solid 
line (NH) and dashed line (IH) distributions. 



C.L. The first set of fits with the values of both CP phases 
assigned to zero, yield the dot-dash line (Amf^ > 0) and 
dotted line (Am^ < 0) distributions. The distributions 
indicate |e eT | < 1.3 at 90% C.L. for either mass hierar- 
chy. For the same scenario but with solar scale mixing 
also neglected (a — > 0), more stringent bounds have been 
obtained by fitting to atmospheric plus K2K neutrino os- 
cillation data [29] . However the bounds become distinctly 
more relaxed when the Dirac CP phase 5 is accounted 
for via marginalization, as shown by solid line (NH) and 
dashed line (IH) distributions from the second set of fits. 
The latter exclusion curves are nearly identical and so 
the hierarchies are not distinguished. At 90% C.L. our 
fits to real-valued e er with S marginalization yield the 
constraint 



- 2.0 < e eT < 2.0 
for either mass hierarchy. 



(29) 



B. 



as sole operative NSI 



For our second scenario we continue to treat e er as the 
sole operative NSI, however we now treat it as a complex 
amplitude by including both its magnitude \e eT \ and its 
CP phase 5 eT in the fit. Additonally we allow the Dirac 
CP phase to be operative. At T2K and MINOS base- 
lines the contribution from solar scale oscillations can be 
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FIG. 2. Allowed-region contours of A^ 2 from fitting to three- 
flavor neutrino oscillations with \e eT \e tS ' T V e as the only op- 
erative NSI matter potential. The regions allowed to |e eT 
and (8 + S eT ) are shown separately for the NH (Fig. 2a) and 
IH (Fig. 2b) neutrino mass hierarchies. In each Figure, val- 
ues within the shaded region bounded by the upper (lower) 
border are allowed by the fit at 90% (68%) C.L. 



regarded as a perturbation, and a simplification arises in 
the limit of no solar scale which can be harnessed to good 
effect. As has been known for some time, the phases S 
and 5 eT , in the limit a — > 0, only appear in the v e appear- 
ance probability as the sum 5 + 8 eT [12] [39] . The verity of 
this assertion can be discerned in part by considering the 
oscillation probability of Eq. ( 26 1 . In Eq. ( 26 ) , the first 
term contains (5 + 5 eT ) as the phase of the amplitude, 



v e ) in which it 
enters~rT 3 T* 2 +T*T 1}2 and |T 3 | 2 ),'the fact that the CP 



/, and the third term contains the same phase combi- 
nation within an oscillatory cosine; otherwise the second 
and fourth terms vanish in the limit a — > 0, and the fifth 
and sixth terms are devoid of phases. 

With careful consideration of the T3 amplitude of 
Eq. (25) and of the terms of Viy^ 

phases only appear in the a — > limit as the sum (8+6 eT ) 
can be seen to hold exactly for v e appearance in constant- 
density matter. 

Since the null solar-scale limit identifies (5 + 5 eT ) to 
be the predominant source of phase in Viy^ — » v e ), we 
express all phases within V(v^ — > v e ) in terms of the sum 
and the difference S±5 eT . We then use the sum-of-phases 
together with \e eT \ as fit parameters, and marginalize over 
the difference-of-phases. As was done for the fit of Fig. [l] 
we also marginalize over the #23 and #13 mixing angles. 
The x 2 fit identifies the regions allowed to the values |e eT | 
and (S + 5 eT ) as shown in Fig. [2] Figure [2^, shows the 
result of fitting to the NH; the result for IH is shown in 
Fig. [2J5. Within each plot, the parameter regions allowed 
by the fit at 68% and 90% C.L. are the shaded areas 
bounded by the lower and upper borders respectively. 

For either hierarchy, there are sizable intervals for the 
sum-of-phases wherein \e eT \ is constrained at 90% C.L. 
to values distinctly smaller than the limit obtained with 
our fit result of Eq. ( 29 ) . The improved constraints for 
\e eT \ in Fig. [2] are made possible by allowing (S + 5 eT ) to 
be a fit parameter; the marginalization of the phase S for 
the fit of Fig. [T] effectively selects phases from regions of 
large excursion in |e eT | as appear in Figs. [2^,, b. 



C. e eT , e ee , e TT with atmospheric constraints 

Given that current limits for e ee and e TT are even less 
stringent than those for e er , full coverage of the possi- 
bilities requires that all three of these NSI be treated as 
operative. Then |e eT |, e ee , e TT , and (S + S eT ) will have sig- 
nificant roles in the fit; on the other hand our data only 
consists of three "bins" of signal rates. We are thus mo- 
tivated to utilize two observations gleaned from analysis 
of this same NSI scenario using the atmospheric neutrino 
data. 

The first observation is that the allowed region of NSI 
couplings is well-characterized by an analytic expres- 
sion [21130]: 



7(i 



(30) 



Relation ( 30 1 is implied by the requirement that oscilla- 
tions with our three NSI couplings be consistent with the 
high-energy atmospheric neutrino data [12) . For our final 
fit we assume relation ( 30 ) to express an equality; with 
this assumption e TT can be expressed in terms of |e eT | and 
e ee , thereby reducing the number of NSI fit parameters. 
Of course, we could as well use Eq. ( 30 ) to eliminate e ee 
instead of e TT ; both approaches have been pursued in the 
literature [rill27] . 
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In the following we use relation (31 ) in conjunction with 
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O Excluded by Atm. v 




exclusion curves obtained from fitting to establish con- 
straints for |e er |. For this purpose relation (31) is very 



useful, for it eliminates a narrow region of otherwise vi- 
able solutions for which (1 + e ee ) ~ and hence e TT , via 
Eq. (30), can be exceedingly large. 

For our final set of x 2 fits, we use \e eT \ and e ee as fit 
parameters. Figure [3] displays our fit results as allowed 
regions in the plane of \e eT \ versus e ee with the distinc- 
tion made between the neutrino mass hierarchies, NH in 
Fig. [3^, versus IH in Fig. |3Jd. The straight-line borders 
of the wedge-shaped region excluded by the atmospheric 
v constraint as encoded by Eq. (31), are superposed on 
each plot of Fig. [3] 

Four separate fits have been carried out in which the 
sum of phases (6 + S eT ) is fixed to a specific value in 
each fit, namely 0, f ,7T, and ~, while the difference in 
CP phases (S — 5 eT ) is marginalized over. The outcomes 
are summarized by the four curves (solid-line, dashed, 
dotted, and dot-dash curves respectively) which appear 
within the shaded areas of each plot. The curves repre- 
sent the boundaries which separate the regions of allowed 
(|e e r|, e ee ) values (areas below the curves) from those 
which are excluded at 90% C.L. These results clearly sug- 
gest that limiting or measuring |e er | at strengths below 
the MSW matter effect is a goal for the longer term. At 
the baselines considered here, even hierarchy discrimina- 
tion in conjunction with a precision S measurement does 
not assure that very restrictive limits are achievable. 

Our most realistic fits however are ones for which both 
of the CP phases are included in the marginalization. 
The outcomes of these latter fits define the parameter 
regions allowed to |e er | and e ee at 68% and 90% C.L. 
as depicted by the shaded areas in Figs. [3}i,b. Thus the 
net effect of the recent T2K and MINOS data is to ex- 
clude those regions of relatively high |e er | which lie above 
the shaded allowed regions and are exterior to the region 
previously disfavored by the atmospheric neutrino data. 



-ee 



VI. DISCUSSION 



FIG. 3. Allowed and excluded regions in the plane of |e eT | ver- 
sus e ee for (a) NH and for (b) IH. Shaded areas denote param- 
eter regions allowed upon fitting to T2K and MINOS data, 
with marginalization over 8 and S eT , as well as over 823 and 
#13. The V-shaped interior region (see Eq. < |31[ )) is excluded 
on the basis of atmospheric neutrino and K2K data [12] . The 
fit sensitivity to CP phases is indicated using separate fits in 
which the value of S + S eT is fixed to, e.g. 0, ^, ir, or These 
give rise to upper boundaries at 90% C.L. for |e eT | as shown 
by the four curves in each plot (solid, dashed, dotted, and 
dot-dashed respectively) . 



The second observation is an approximate bound de- 
duced from atmospheric neutrino data [12\ 130) : 

I CeT I < 1.1 X |1 + e ee \. (31) 



Previous investigations of NSI matter effects in neu- 
trino oscillations were hindered by lack of a measured 
value for the #13 mixing angle. The present work has 
availed itself of the recent delineation of #13; it is the 
first study to use v e and V e appearance measurements 
from accelerator long-baseline experiments, in conjunc- 
tion with atmospheric-neutrino measurements, to obtain 
constraints for the complex e er NSI coupling. The con- 
straints are expressed by the allowed regions in Figures [2] 
and [3] The limiting values for | e er | vary according to the 
sum of the CP-violating phases 6 + S eT and according to 
the choice of neutrino mass hierarchy. At 90% C.L. the 
maximum value allowed to |e eT | varies from 0.7 to 2.3, 
as shown by the minima and maxima in the allowed re- 
gions versus 6 + S eT for each mass hierarchy. These values 
represent an improvement upon the limit \e eT \ < 3.0 pre- 



viously inferred from world data [14] . Nevertheless, the 
allowed range for \e eT \ permitted by our analysis is still 
relatively large, with the strength of the NSI potential 
exceeding that of the conventional MSW matter effect 
remaining as a viable possibility. Our fit results of Fig. [3] 
show that, with the flavor-diagonal NSI e TT expressed 
in terms of |e eT | and e ee using the atmospheric relation 



(301, the v e appearance data does not provide any upper 
limits to the magnitude of the e ee NSI. A v e appearance 
measurement to constrain e ee will likely require signifi- 
cantly more data, obtained with detectors located at two 
different baselines such as in the 295 km and 1050 km 
baselines of the T2KK proposal [H] . 

In the near term, the T2K experiment is pursuing a 



precision measurement of Viy^ 



with Vn expo- 



sures continuing at higher beam power. A v e appear- 
ance sample with sevenfold more events is projected for 
2015 gl] ■ For the MINOS experiment on the other hand, 
running with the NuMI neutrino beam in its low-energy 
configuration has completed. A new round of data tak- 
ing with the MINOS detectors will commence in 2013 
with the NuMI beam operating in medium-energy mode, 
for the MINOS+ experiment. With medium-energy run- 
ning, neutral-current interactions yielding shower-like fi- 
nal states will occur at higher rates than was the case in 
MINOS exposures and so v e appearance measurements in 
MINOS+ may not be feasible. Medium-energy running 
of the NuMI beam however is optimal for v e appearance 
measurements using NO^A. The 810-km long baseline 
of NOf A is off-axis with respect to the NuMI beam; 
it receives a narrow-band neutrino flux for which back- 
grounds originating from high energy neutral-current re- 
actions are mostly suppressed. It is the NO^A experi- 
ment whose observations of v e /v e appearance rates and 
spectra, when taken together with new T2K measure- 
ments, hold promise for significant gains in delimiting 
NSI matter effects for neutrinos in propagation [T7] . The 
discovery reach for NSI of the projected T2K and NO^A 
exposures, when combined with measurements from the 
reactor experiments, has been examined in Ref. [30] • A 
point often made is that a degree of redundancy among 
the international suite of neutrino long-baseline acceler- 
ator and reactor experiments is useful for delineating the 
standard three-flavor oscillation framework. With NSI 
matter effects for neutrinos in propagation included for 
consideration, multiple measurements conducted at dif- 
ferent baselines will be essential to affirming or ruling out 
the e eT non-standard interaction. 
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APPENDIX: DERIVATION OF Afa -> v e ) 
INCLUDING e er , e ee , AND e TT NSI 

Our derivation of the exact A(v^ — > v e ) for the matter 
Hamiltonian of Eq. pi for neutrinos propagating through 
a constant-density medium, proceeds as described in Sec. 
IV of Ref. [31] . In brief, the conventional three-flavor 
Hamiltonian in flavor basis is transformed to the 

propagation basis H ^ . Upon re-phasing of the diagonal 
elements of (with minor differences from the de- 

scription in Ref. [31]). one arrives at the Hamiltonian of 
Eq. (13). We separate into an "unperturbed" part, 



Hq , plus an interaction potential, 



(p) 



-Q 









r 







-G 




h 







/* 





+Q 




b* 


: 



H, 



We then define an Interaction Picture: 



(32) 



j7«(t) = e iif o V p) (*), v (p) (t) = e~ iH ° WHt) , (33) 



so that 



where 



i — v 
dt 



(t) = V I -^ I \t) 



Vi(t) 



iff 



-iff 



(34) 



(35) 



Our approach is to solve for the time evolution opera- 
tor in the Interaction Picture Ui(t = £, t = 0): 



^ I \t) = U I {t,Q)-0 {I \O). 



(36) 



Substitution of Eq. (36) into Eq. (34) yields the wave 



equation which governs Ui(t, 0) 
. d 



%-U I (t,0) = V I (t)U I (t,0). 
at 



(37) 



To obtain Vi(t) we require the matrix representation 
(in propagation basis) of the unitary operator forms 

exp(±iH^ p) t). From H^ p) we extract the reduced matrix 



r>(p) _ 
H o,R ~ 



-Q fo + Hi 
fo ~ ifi +Q 



(38) 



where /o and j\ designate the real and imaginary parts 
of the element (i?Q P ')i3 = /. The reduced matrix can be 
decomposed using Pauli matrices, 



H, 



(p) 



0,7? ~ fo^ x ~ f^V ~ Q® z 

= N-a, where N = (/ , -/i, -Q). 



(39) 



We have |jV| = iV = y/f§ + fl + Q 2 . Its unit vector h 
defines the axis of rotation in the reduced (spinor) space, 



1 



tt(/o,-/i,-Q). (40) 
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Designating the angle of rotation with 

(j> = Nt, 
we use the spinor identity 

ia-n<f> _ ( cos (f> + in z sin (f) (in x + n y ) • sin ( 



(in x — riy) ■ sin <f> cos <j> — in z sin <^> 

and furthermore define 
7 = cos + in z sin <^>, 

,S = P x + i(3 y , P x = n x sin cj>, j3 y = n y sin <fi. 



(41) 



(42) 



(43) 



Then i/3 = i/^ — P y and = + P y , and we have 



7 iP* 
ip 7 * 



Thus in the propagation basis we may write 



e o 



7 if3* 
e- iGl 
ip 7* 



(45) 



To move the formalism to the Interaction Picture, we 
evaluate 



Vi{Z) = e lH ° *• V-e-'" 1 

u 
= | u* u 

v* 



where the complex elements of Vi(£) are 

u = (tt + i(3*B)e tGe , v = foB - ip*r*)e~ iGt . 
Now 



(46) 



|u| 
|u| 2 + |w| 2 



(47) 



(48) 



(m;)* 







The real- valued expression (|u| 2 + |i>| 2 ) recurs upon taking 
higher integer powers of Vi(£). It is readily reduced to 
(|r| 2 + |&| 2 ), previously designated as rj 2 in Eq. (21): 



rf = \u\ 2 + \v\ 2 = \r\ 2 + \bf 



(49) 



The exponentiation of Vi(£) into e proceeds as in 
Ref. EH. We obtain 



We define 



Vi 



(1 - cos(r)£)) - i— sm(r)£). (50) 
rj 1 rj 



- - U r- V 
T)£, u = -, V= -, 



(51) 



and write the diagonal elements of Eq. ( 50 1 as 



D u = l-2|u| 2 - sin 2 d=cos6, D v = 1 - 2|v| 2 • sin 2 -. 

2 (52) 
For the off-diagonal elements we define 



w — u sin 9, p—~2uvs'm — , k = vsin( 



(53) 



Then the evolution operator in the Interaction Picture is 

(54) 



D u -iw p 
Ui(£,Q) = e~ lVllL = | -iw* d -i 
p* -ik* D, 



and, in the propagation basis, it becomes 



(44) U { p\£,0) = e- i6 ° )t -U I (£,0) 



-y*D u -iP*p* 7*(-w) - P*k* ~f*p-i/3*D v 



-iw e 



iGf 



de 



iGl 



(—ik)e 



iGC 



7J3* — i(3D u "/(—ik*) — f3w "/D v — i(3p 



(55) 



Finally, returning to flavor basis 

(£, 0) = ■ (£, 0) • (LsRf) (56) 

we obtain 

C23C / i2 + S23(J 13 e — 

{-i)c2z{l* w ~ ifi*k*) + s 2 z~l*pe~ lS - is 2 zP* D v e~ lS . 

(57) 



We insert expressions ( 52 ) and ( 53 1 for the elements of 



k, p, and D v into the last line of Eq. (57 1 and re-arrange 
the order of the terms to obtain 

My* -> v e ) = H)s2 3 P*e- lS 



+ (~i)c23 [l*u- i(3*v*]sm6 

+ 2s 23 [i/3*|u| 2 - l*uv] • sin 2 - • e - *. 



(58) 



The three terms of Eq. ( 58 ) correspond, respectively, 
to the terms T 1 +T 2 + T 3 of Eq. |(9j). 



For the first term, T l7 we use Eqs. (43), (41 ), and (40) 
to write 

p* = (n x - in y ) sin(iW) = i(/ + i/i) sin(^); (59) 



the term reduces immediately to 



/ 



T 1 = H) S23 ^sin(iVA)- e - 



(60) 



Considering the second term T 2 , we use Eqs. (147b and 



(51 1 to insert 

u = -for + iP*B)e lG \ v* = -( 7 *B + i/3r)e lGf , 
»7 ^7 
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and find that it reduces to 



Substitution of 7 and /3* from Eq. ( 43 ) yields 



T 2 = H)c 23 - -sin^ A) -e l0A . 
V 



(61) 



The remaining term is 



T 3 = 2s 23 [if3*\v\ 2 - 7 *uv] • sin 2 - • e~ l8 . (62) 



The expression within the bracket reduces to 
4 [--/rb + i(3*\r\ 2 ] so that 



(63) 



T3 = (-2*23) • sin 2 - • e 



rb 



\r\ 2 n y 
i} 2 



rb 



sin <p 
(64) 



Within the curly brackets on the right-hand side, the 
factors multiplying cos <j> comprise the complex function 
Si of Eq. (23), and the expression which multiplies sin (9 



is the complex function 5 2 of Eq. ( 24 ) . Thus Eq. ( 64 ) 
coincides with Eq. (25) for T3. 



With Eqs. (60), (61), and (64) we have shown that 
Xi, T 2 , and T3 have the forms as previously specified 
in Eqs. (pi), pOL and p5|) of Sec. ImBl The transition 



amplitude Aiy^ 
ified. 



v e ) of Eq. ([9| is thus completely spec- 
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